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ABSTRACT 



While the average metallicity of the intergalactic medium rises above Z > 
1Q~^Zq by the end of the reionization, pockets of metal-free gas can still exist at 
later times. We quantify the presence of a long tail in the formation rate of metal- 
free halos during late stages of reionization (redshift z ^ 6), which might offer 
the best window to detect Population III stars. Using cosmological simulations 
for the growth of dark matter halos, coupled with analytical recipes for the metal 
enrichment of their interstellar medium, we show that pockets of metal-free gas 
exist at ;z ~ 6 even under the assumption of high efficiency in metal pollution via 
winds. A comoving metal-free halo formation rate d'^n/dtdV > 10^^ Mpc^'^ yr^^ 
is expected at 2 = 6 for halos with virial temperature T^r ~ 10^ K (mass 
~ 10^ Mq), sufficient to initiate cooling even with strong negative radiative 
feedback. These halos will appear as absorption systems with a typical hydrogen 
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column density of ~ 10^° cm~^, a sky covering factor 5 x 10~^ and a number 
density of 25 arcmin~^ for 5.5 < z < 6.5. Under the assumption of a single 
Population III supernova formed per metal-free halo, we expect an observed 
supernova rate of 2.6 x 10"'^ deg~^ yr~^ in the same redshift range. These metal- 
free stars and their supernovae will be isolated and outside galaxies (at distances 
> 150 kpc) and thus significantly less biased than the general population of 
~ 10^ M0 halos at 2; ~ 6. Supernova searches for metal-free explosions must thus 
rely on large area surveys. If metal-free stars produce very luminous supernovae, 
like SN2006gy, then a multi-epoch survey reaching mAB = 27 at 1 /im is sufficient 
for detecting them at 2; = 6. While the Large Synoptic Survey Telescope (LSST) 
will not reach this depth in the z band, it will be able to detect several tens of 
Population III supernovae in the i and r bands at z < 5.5, when their observed 
rate is down to 3 — 8 x 10~^ deg~^ yr~^. 

Subject headings: galaxies: high-redshift - early universe - intergalactic medium 
- supernovae: general - cosmology: theory - stars: formation 



Introduction 



Metal-free stars, conventionally called Population III stars, represent the first step to 
exit the dark ages of the universe and start the reionization era. The chemical enrichment 
of the intergalactic medium (IGM) opens the way t o the formation of stars and galaxies 
like our own Milky Way (see iBromm fc Larson! |2004| ). This first step in the star formation 
history of the Universe has attracted large theoretical interest, and thanks to a combination 
of analytic and numerical investigations, we now have a reasonably solid understanding of 
the properties of Population III stars within the ACDM paradigm. Population III stars 
are expected to start forming within dark ma tter halos of mass ~ 10^ M^r^ (minihalos ) 



at z ^ 60, about 40 Myr after the Big-Bang jNaoz et all bood : iTrenti fc Stiavelh 



most likely with a top-heavy mass function and just one star per halo (lAbel et al. 



2007h 



2002 



Bromm fc Larson I2OO4I : lO'Shea fc Norman! bopj boosh . Once the epoch of Population III 
stars begins, their formation rate in minihalos increases exponentially with redshift, reaching 
a level of 10"'' Mpc~'^ yr~^ at 2; ^ 25. Eventually, by 2; ^ 15, Population III star formation in 
minihalos is suppressed by the build-up of a strong radiative background in the II2 Lyman- 
Werner (LW) bands, which photodissociates the molecular hydrogen responsible for the 
cooling of the primordial gas. A second mode of metal-free star formation is possible in 
larger halos (M > 2 x 10'' M0) where self-shielding against LW feedback is more effective 
and where cooling may also be driven by Lya lines. However as the redshift decreases, it 
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beco mes progressively more un l ikely that ~ 10^ Mf;^ ha los composed of metal-free gas can 
form (IFurlanetto fc Loebl l2005l : iTrenti fc Stiavellil |2009| ) . The transit ion to Population II 
star formation happens in the metallicit y range Z ~ 10~^ — 10~^ Zq (IBromm et al.ll2001al : 
Santoro fc ShuU 20061: Smith et al. 2009 ) or even at lower metallicity in the presence of dust 
( ISchneider et al.ll2006l ). Pockets of pristine gas are still likely to exist at z < 6, thanks to the 



inhomogeneity of metal en r ichment in the IG M (jScannapieco et al.l l2003l : ISchneider et al. 
20061 : iTornatore et~aDl2007l : iRicotti et allboosh . 



An unambiguous observational test of the Population III paradigm would represent 
a significant step toward our understanding of star formation in the Universe. This is 
still missing, despite many attempts ranging from Galactic su rveys for metal-poor stars 
( iBeers fc Christliebl l2005l : iFrebel et al.l 120071 : ICohen et al.l l2008l ) t o surveys aimed a t iden - 
tifying their spectral signatures in galaxies at very high redshift ( iNagao et al.ll2005l . l2008l ). 
In fact, Populat ion III stars are expected to be extreme ly hot and thus strong emitters in 



the He II lines ( iTumlinson fc ShuU 



200G 



Oh et a 



been proposed. For example, iJimenez fc HaimanI (120061 ) argue that many puzzles related 



2001 ). Some indirect evidence has also 



to Lyman Break Galaxies at 2; > 4 would be solved by assuming t hat a fraction of their 
light is produced by metal-free stars. Similarly, iMalhotra fc Rhoadsl (120021 ) suggest that the 
l arge equivalent width of Lya emitters ed. z > 5 might be a signature of Population III stars 
( ITumlinson et al.ll2003l : IScannapieco et al.ll2003l ). However co-existence of both Population 
HI and II stars in the same galaxy appears to require a fine tuning of the star formation 
history, especially because the initial mass function (IMF) of metal-free stars is probably 
more biased toward large masses than that of Population II stars. In addition, it would 
be extremely hard to avoid mixing of of metal-fr ee and metal polluted gas in a galaxy, be - 



cause of the turbolent nature of the ISM (e.g. see lBalsara fc Kimll2005l : iBalsara et al.ll2008l ). 
Therefore a significant population of primordial stars in Lyman Break Galaxies is unlikely. 
Another indirect fingerprint of Population HI stars is the presence of peculiar abundance 
patterns in the gas cher nically enriched by Population HI stars, ^vith en hanced even-odd ele- 
ment abundance ratios (jHeger &: Woosleyll2002l : ITumlinson et al.ll2004j ). but again no strong 



evidence is available. Finally, the transition from Population HI t o Population II star forma- 



tion lea ves an imprint in t he reionization history of the Universe (jShuU fc VenkatesanI 12008 
see also lBagla et al.ll2009l ). 



The null result achieved so far is not surprising. In fact, the task is very challenging: 
the expected top-heavy IMF implies that the lifetime of a single Population HI star is 
rather short — possibly only a few milli on years — and they form in isolation, at least in 



minihalos (e.g., lO'Shea fc Normanll2008l ). These objects are extremely faint at cosmological 



distances (mAB > 40 at 2 > 10) and thus unlikely to be detected directly via imaging even 
with future facilities such as the James Webb Space Telescope (JWST). A better chance 
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of detection might be afforded by pair-instability supernovae. In fact, metal-free stars in 
the mass range 140 — 260 Mq end their lives with powerful explosions (up to 10^^ ergs 



released) powered by the production of e~^e~ pairs (IHeger &: Woosleyl |2002| ) . These are 
up to 10 times more lum inous than standard SNs from Population III stars of lower mass 
( iHeger fc Woosleyl |2002| ). although also metal-free stars with ~ 30 — 100 Mp-, might p roduce 
luminous hypernovae with an energy release above 10^^ ergs (INomoto et al.l 120031 ). Still, 
many observational challenges are present at such high red shift (z > 10), especially because 
of the significant tirn e dilation of their light curves (e.g. see IScannapieco et al.ll2005l : see also 
Mesinger et aPbooeh . 



In this paper, we aim to evaluate the number density of pockets of primordial gas that 
survive to z < 10, using conservative assumptions on the metal enrichment efficiency and 
parameterizing the results based on the strength of galactic winds. Even with efficient metal 
transport, metal-free pockets still survive at z ^ 6. Thus if massive stars are formed out of 
this primordial gas, the expected star formation rate can lead to the detection of hypernovae 
using all-sky surveys such as the Large Synoptic Survey Telescope (LSST). 

The paper is organized as follows. In §[2] we present our model for the formation rate of 
metal-free halos around the end of the reionization era, whose predictions are discussed in 
§ |3] along with the prospects for the observational measures of the metal-free star formation 
rate after reionization. We conclude in § |H 



Pristine halo formation rate 



2.1. Numerical setup 



The numerical simulations have been carried out using the Particle Mesh Tree code 
Gadget- 2 (Springei 2005 ). We adopt a cosmology based on the fifth year WMAP data 



komatsu et aPboOsh 



20081 1: Qa = 0.72, Qrn = 0.28, Qb = 0.0462, Hq = 70 km s'^ Mpc"^ and 
(Js = 0.817, where flm is the total matter density and fib is the baryon matter density in units 
of the critical density, pc = 3Hq/{87iG). Here, Hq is the Hubble constant parameterized 
as Hq = 100 h km s~^ Mpc~^, Q\ is the dark energy density, and Ug is the root mean 
squared mass fluctuation in a sphere of radius 8 h^^ Mpc extrapolated to 2; = using 



algorithm ( 


Bertschineer 


Eisenstein & Hu 


( 


199d) 



0.96. 



N 



We start our principal simulation at redshift z = 199 using a box of edge 1 = 7 h-^ Mpc, 
= 1024^ particles, a mass resolution of 3.4 x 10^ M© and a force resolution of 0.16 h^^ kpc. 
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A total of 79 snapshots are saved from z = 30 to z = 4, and from these snapsh ots dark 
matter halos are identified with a friend-of-friend halo finder (jPavis et al.l Il985l ) using a 
linking length equal to 0.2 of the mean particle separation. An additional, lower resolution 
simulation (A^ = 512^, / = 10 Mpc) is described in Appendix [XJ where we discuss the 
convergence of our results. 

For every redshift z, we define a minimum halo mass for primordial star formation which 
is the lowest halo mass that can efficiently cool in a Hubble time. Two ma in coolants are 



available for metal-free gas: molecular hydrogen and atomic hydrogen (see iTegmark et al. 



19971 ) with different efficiencies depending on the gas temperature. At the time of formation 
of the halo the gas temperature is approximately equal to the halo virial temperature, 



Tvir(M,z) = 2554 K 



M 



106 Mq 



2/3 



1 + Z 

31 



M olecular hydrogen is an effective coolant at virial temperatures Tyir ^ 10^ K ( Galli &: Palla 
but it can be photodissociated by radiation in the H2 Lyman- Werner bands (iLepp fc Shulll 
1983. ) ■ which progressively builds up as the redshift decreases and the Universe starts to be 
ionized. Extended P ress-Schechter modeling shows that this cooling channel is significa ntly 
quenched at 2; > 15 (IHaiman et al.lll997l : iHaiman fc Bryanll2006l : iTrenti fc Stiavellill2009l ). at 
which point metal-free star formation shifts to halos with T^i,. > 10^ K, where self-shielding 
from the Lyman- Werner background is efficient and atomic hydrogen cooling is also possible. 
We therefore assun ae that the minimuna hal o mass capable of cooling is that derived in the 
reference model of iTrenti &: Stiavellil (120091 ) . which includes the effect of a self-consistent 
Lyman- Werner background. This minimum mass versus redshift, M^i^{z), is presented in 
Fig. [H Given the mass resolution of our simulation, we can resolve halos capable of forming 
Population III stars via H2 cooling up to 2; < 26 with A^ > 32 particles. For 2; < 13 the 
minimum halo mass capable of forming stars corresponds to that of halos with Tvir > 10^ K, 
which can be written as: 



Mrniniz < 13) = (7.75 X lO*' Mq) 



1 + Z 

31 



-3/2 



(2) 



As a conservative upper limit to the self-enrichment of halos, we flag all particles of a 
halo with mass M > M^ij^{z) as "metal enriched". We further assume that a necessary, but 
not sufficient condition for a halo to be pristine is to be made entirely of particles that have 
never been in a halo capable of cooling at higher redshift. 

Metal enrichment is also possible due to supernovae winds originating in nearby halos. 
Wind bubbles during the reionization epoch can propagate out to a comoving radius of 
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~ 150 kpc by redshift z = 6 accordingly to the iFurlanetto &: Loebl (120031 ) model (but 
Madau et al.ll200ll predict bubble radii approximately two times smaller). An estimate of 
the wind bubble radius ca n be also obtained noting that a bubble expanding at ~ 60 km s~^ 
(IFurlanetto fc Loebl l2003l ) reaches a radius of ~ 150 kpc after 500 Myr. Metals are in 
fact assumed to be transported by winds at high speed (~ 300 km s~^) until they mix in 
a thin spherically symmetric shell, propagating at ~ 60 kn i s~^, that eventually fragments 
by cooling and pollutes the IGM (IFurlanetto fc Loebl l2003l ) . The simple picture discussed 
here is consistent with more detailed modeling of the interaction of the win ds with the IGM. 



Using the non-cosmological Sedov- Taylor model, as written in equation 8 of iTumlinson et al. 
(120041 ). we obtain bubble sizes of < 150 kpc at z ~ 6. 



To derive a conservative upper limit to the wind enrichment, we compute for every halo 
that is not self-enriched the distance to neighbor halos that have experienced star formation 
in the past, that is of mass M > Mmin. The size of the wind bubble around these halos 
depends on the first time they experienced star formation, assuming a bubble expansion 
speed in the range fwind = 40 — 100 km s~^. If the distance between the non-self-enriched 
halo and other halos is greater than the wind bubble sizes, we have a truly pristine halo. No 
m etal wind had the op portunity to reach the halo (see also the discussion of metal enrichment 



m 



Ricotti et al. 20081 



This assumption is in reality likely to overestimate the efficiency of wind pollution. As 
the wind bubble expands, its capability of enriching the IGM above the critical metallicity 
needed to quench Population III formation is progressively reduced. For example, if we 
assume a critical metallicity of 10~^'^ Zq, a non self-enriched halo of mass 10^ Mq needs to 
be polluted with ~ 100 Mq of metals by winds. We assume that the galaxy with the winds 
has a star formation efficiency e = 10~^ and a metal yield of F ~ 2.4 x 10^^ (total mass 
of metals released per unit solar mass in the IMF). This value is estimated by assuming a 
3 Mq metal release from a supernova and one supernova formed for every 120 Mq of mass 
in the IMF. The wind then contains about 3960 Mq of metals for a wind galaxy also hosted 
in a 10® Mq halo. Therefore a non self-enriched halo needs to have a cross section with the 
wind bubble of at least 100/3960 x An ^ 0.32 sr in order to be polluted above the critical 
metallicity assuming a 100% efficiency in the pollution process. Given that a 10® Mq halo has 
a (comoving) virial radius ~ 10.5 kpc at z = 6, the maximum radius of the wind bubble 
is only ~ 33 h^^ kpc. Of course, multiple bubbles can contribute to the metal enrichment of 
a halo, and this enrichment can happen before the virialization process is completed. The 
gas is then less concentrated and has a larger cross section. 
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2.2. Results 



In Fig. [2] we plot the formation rate {(Pn/dtdV) of halos that are not self-enriched. 
This quantity is slowly declining from z = 25 to z = 13 and exhibits a local maximum at 
z ~ 10, at which time about 30 halos are formed per year per comoving Gpc^. Interestingly, 
their formation rate is only slowly varying in the redshift range 7 < 2; < 25 because of 
two competing effects. First, overdensities that lead to these halos become progressively 
more efficient, and second, as time passes, the supply of metal-free gas decreases and it is 
more probable that a region has been self-enriched. The halo formation rate starts dropping 
rapidly only at 2; < 7. The evol ution of cP"nldtdV measured in our numerical simulation 
is in excellent agreement with the lTrenti fc Stiavellil (120091 ) model, shown as the red-dashed 
line in Fig. [2j There is a discrepancy by a factor ~ 2 at 2; > 15, probably related to the finite 
resolution of the numerical simulation which does not allow us to resolve the progenitor of 
halos formed at these redshifts. The agreement improves to ~ 30% at lower redshift. With 
the aid of the analytic model, we can u nderstand the origin of t he local maximum in the halo 
formation rate at 2; ~ 10 (see Fig. 1 in lTrenti fc Stiavellil l2009l ) as associated with the peak 
of Population III stars formed in halos with Tvir = 10^ K. The decline from 2; ~ 25 to ^ ~ 13 
is instead associated with Population III formation in minihalos, progressively quenched by 
the rising Lyman- Werner background. 

The spatial distribution of these halos is consistent with that of the general population 
of halos with the same mass, independent of their enrichment history, as we have verified by 
measuring the two-point correlation function and thus their bias. This result is natural in 
the context of extended Press-Schechter modeling, because self-enrichment depends only on 
density fiuctuations at scales smaller than the halo mass. 

Dependence on the environment is present instead in the case of wind enrichment. To 
illustrate the effect of wind enrichment, we plot in Fig. [3] a distance vs. mass scatter plot for 
the nearest neighbors of non self-enriched halos formed at 2; = 6. The halos closest to those 
that have not been self-enriched are typically small (10^ M©) and at distances comparable 
to those over which the wind bubbles can propagate at 2; = 6. Metal enrichment of non-self 



(Stocke et al. 


2006 


; see also 


Fujita et al. 


— — 
2004) 



halos have no neighbors within 150 h^^ kpc, a distance that essentially guarantees that these 
halos contain genuine metal-free gas. Depending on details of the wind propagation, up to 
50% of these halos can be intrinsically metal free. The median of the distance distribution 
in Fig. [3] is 82.4 kpc, for which a single wind bubble is unlikely to chemically enrich 
a metal-free halo above the critical metallicity for quenching Population III formation. To 
quantify the wind enrichment we compute the wind bubble size for all the star-forming halos 
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in our simulation as described in § 12. 1[ We plot the resulting metal-free halo formation rate 
in Fig. Hlfor three different values of the wind bubble expansion velocity, fwind = 40, 60, 100 
km s~^. Even in the "worst case scenario" of propagation at 100 km s~^, metal-free halos 
continue to exist down to 2; ~ 5. 



2.3. Absorption Systems Predictions 



These metal-free halos leave an observational signature as absorption systems before star 
formation is associated with them. To estimate their column density, we assume that the 
baryon profile follows the dark matter, a reasonable assumption immediately after virializa- 
tion and before significant c ooling takes place. A ssuming a helium mass fraction Y = 0.2477 
and a WMAP-5 cosmology (IKomatsu et al.ll2008l ). we can write the typical hydrogen column 
density of an halo of total mass Mh as: 



(Nh) 



10^° cm-2 



n. 



0.1372/ VIO8 Mo 



1/3 



l + z 
7 



(3) 



If the gas spends 5 Myr before cooling and initiating the protostellar core formation, we 



estimate from Fig. H] that the number density of metal-free absorbers is 
z = 6. Their physical size is 

1/3 / 



3M, 



h 



4vrp^ 



;i5.3 kpc)(l + 2) 



h 



108 Mq 



25 arcmin ^ at 



(4) 



and their covering factor on the sky is ~ 5 x 10"'^ for z G [5.5 : 6.5]. The column density 
of these system is at the edge of the Damped Lya (DLA) limit and thus well within the 
observational limits of both Lya absorption or 21-cm emission studies. The main limitation 
for the detection in absorption is the very low number of bright QSOs known at 2; > 6 
( iFan et al.l l2006l ). In addition, the number density of metal-enriched DLAs is higher by 
a factor of ~ 300 compared to that of metal-free systems if we extrapolate to z ~ 6 the 
intrinsic density of neu tral gas found at z ~ 4.5 in DLAs {n^Hi^^\z = 4.5) ~ 10-3; see 
Trenti &: Stiavellil l2006l ) Obtaining strong limits on the metallicity of a z ~ 6 absorber will 
be challenging even for a 30 m class telescope. 



3. Population III stars and prospects for their detectability 



From the metal-free halo formation rate we can derive a Population III star formation 
rate once we assume an efficiency for the conversion of gas in stars, a poorly constrained 
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parameter. Numerical simulation of the collapse and cooling of metal-free halos seem to 
favor massive Population III stars formed in isolation, even in halos with virial temperature 
T yjr ~ 10^ K in the presence of a strong Lyman-Werner background (see lO'Shea &: Norman 



20081 ). However, these simulations are limited to the formation of protostellar cores and 



are not able to reach the point of actual creation of the star. In addition, cosmological 
simulations focus on Population III star formation at the start of reionization at z > 15. 
The physical conditions at the end of the reionization epoch are different. In particular, the 
CMB temperature floor is significantly lower, and a strong ionizing UV background is present, 
wh ich may lead t o a cha nge in the initial mass function of Population III stars as suggested 
by lYoshida et al.l (120071 ). Late time metal-free stars could have masses of 30-50 Mq, rather 
than the values M > lOOM© suggested at higher redshift. Note also that, while the very 
massive star hypothesis is a strong prediction of numerical simulations of Population III star 
formation, there is some tension with the mass function neede d to explain the abundance 
patterns of extremely metal poor stars in the Milky Way (see iTumlinsoru 120061 ). Here we 
explore these two different scenarios for the formation of metal-free stars at the edge of 
the reionization, by assuming (i) very massive stars formed in isolation and (ii) clusters of 
massive stars. 



3.1. Detection during the main sequence lifetime 

First, we evaluate the luminosity of Population HI stars on the main sequence. Under 
the assumption that these s tars are massiye, we find their main sequence luminosity to be 



3 X 10^^ erg s ^ Hz ^ Mq^ ( Bromm et al.l l2001bl ). which translates into a 1 /xm observed 



magnitude ~ 40 at 2; = 6 for a 100 Mq star. Therefore, direct detection will be 

extremely challenging (see Tabled]). This is essentially hopeless if star formation happens 
in isolation. Even if the massive star formation formation occurs in clusters, significant 
problems are present. For example, a cluster formed with star formation efficiency e = 10~^ 
in a 10^ Mq halo a would have an observed 1 fim luminosity ~ 32. 

In our model with fwind = 60 km s^^, we expect a number density on the sky of ~ 10 
arcmin"^ i-dropout^ Population HI sources at 5.5 ^ z < 6.5 (assuming that each source 
remains visible for 2 Myr). Detecting them requires significantly deeper exposures than th e 



HST-UDF, as those observations only reach down to ^ssolp = 28.6 at 5a (lOesch et al.ll2009l ). 



Even JWST might not be sufficient, unless the Population HI star formation efficiency is 



^ These are high- redshif t sources identified in broad-band imaging by inferring the location of the Lya 
IGM-absorption break, see ISteidel et al.l (|l996l ). 
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higher than e = 10~^ (see also I Johnson et al.ll2009l ). In fact, a deep-field JWST integration 
is expected to reach ^ 31 in 65 hours (see Table [1]). Adopting a different strategy 
and searching for Lya emission vi a narrow-band filters would require a sensitivity of 3 x 
2q4i gj,g g-i ^Qj, ^Y\_e same sources ( Bromm et al. 2001bl ) with an expected number density 
~ 0.5 arcmin"^. The lower number density on the sky is given by the smaller redshift range 
probed by narrow band filters {Az ~ 0.05 versus Az ~ 1 in broad-band imaging). Current 8 
m class telescopes reach a s ensitivity of 10^ ^ erg s~^ at 2; ~ 7, equivalent to an observed fiux 



of 2 X 10 erg s ^ cm ^ (lOta et al.l 120081 ). A 30 m class telescope is expected to improve 



this limit by about a factor 50, reaching ~ 10 erg s 1 cm , which might be sufficient to 
detect Lya emission from a large cluster of Population III stars. 

Even if future telescopes reach the sensitivity required to detect Population III clusters, 
the main problem is how to distinguish them from Population II so urces. Adopting the 
luminosity function measured at redshift z = 6 from the UDF survey (IBouwens et al.ll2006l ) 



and extrapolating it down to = 32, we expect a number density of ~ 42 arcmin"^ i- 
dropouts from metal-enriched galaxies. While contamination from Population II stars can 
i n principle be addressed owing to the higher effective temperature of Population III stars 
(ITumlinson fc Shulll I2OOOI ) . it will be extremely challenging to identify the characteristic 
He II emission features expected from higher temperatures, for sources near the detection 
limit of an imaging-only survey. Therefore, disentangling the relative contribution does not 
appear feasible, especially given the int rinsic uncertainties on the luminosity function at 
high redshift (see lTrenti fc Stiavellill2008l ). Although Population III sources at lower redshift 
[z ~ 5) will appear brighter than their z ^ 6 counterparts, their direct identification as 
V-dropouts will be even more challenging because of the increased number of Population II 
V-dropout sources combined with a sharp decline of the number of Population III sources 
at this redshift (see Fig. H]). 



3.2. Supernovae Detection 



The possibility of detecting Population III stars in their Supernova phase are more 
promising, especially if these stars are very massive and hypernovae and/or pair instability 
SNe (PISNe) are prefere ntially produced. Theoretical modeling of PISNe light curves by 
Scannapieco et al.l (120051 ) predicts a peak luminosity Mab G [—19 : —21.5], depending on 
the mass of the progenitor (M G [140 : 26O]M0). Thi s is in reasonably g ood agreement with 
the most luminou s SN o bser ved to date, SN200 6 gv ([Smith et al.l 120071 ). but the results by 
Heger fc WooslevI J2002h and IScannapieco et al.l J2005f i" are based on one-dimensional, non- 
rotating simulations. Thus they carry some uncertainty both in the expected mass range 
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for PI SN expfosions as well as in the total energy released. For example, lEkstrom et al. 



( 120081 ) find that inclusion of significant amounts of rotation could enhance the explosion 
energy. Assuming here the observed light curve of SN2006gy as a reference (but note that 
SN2006gy is not metal free), we expect a peak luminosity Mab ~ —21.8, which translates 
into a 0.14 fim Mar ~ —20.5 assuming a blackbody spectrum with effective temperature 



Teff = 15000 K (jScannapieco et al.ll2005l ). This leads to an observed 1 /xm mAB ~ 26.25 at 
z ~ 6 (distance modulus of 48.86 and K-correction of — 2.51og^o (1 + ^))- To obtain a robust 
detection of the supernova peak luminosity curve, a survey reaching down to rriAB ~ 27.0 
might be adequate, although the near-UV spectrum of SN2006gy is not well constrained. 
Thus the prediction of its observed luminosity is intrinsically uncertain at A < 1/im when 
the source is at z ~ 6. Because a pair instability supernova light curve remains near its 
peak for a relatively long time (two to three months) plus cosmological time dilation for the 
observer rest frame, a search for z ~ 6 SNs would require a multi-year time-frame. This is 
convenient, as it allows stacking multiple i mages for recurrin g all-sky surveys such as LSST. 



Given the current LSST requirements (see llvezic et al.ll2008l ) a 5a detection at ttiab = 26.25 



in the z band can not be reached wi thin the main survey, but is within reach of a deeper 



■ mim-survey (see 



Ivezic et al.ll2008l Sec. 3.2.2). 



Supernovae detection at 2; ~ 5 does appear potentially within the capabilities of the 
main LSST survey, even though these sources will be rarer than their z ~ 6 counterparts. 
Population III supernovae at 4.5 < < 5.5 can be detected with LSST in the i-band, where 
the instrument is more sensitive by 0.7 magnitudes than in the z-band. In addition, the 
z ~ 5 sources will be about 0.3 magnitudes brighter compared to those at z ~ 6. A Aa limit 
of rrii = 25.95 can be achieved by stacking the 23 individual exposures, each with sensitivity 
of 24.0, that are planned for each year of operations (see Table [1]). Combining the i-band 
data with those in the z-band will further increase the confidence level at which a source 
is identified. This strategy will push the observations to their detection limit, but if there 
is a distribution of peak luminosity, LSST may be able to identify at least the brightest 
explosions. Of course, with the very large-data data set generated by LSST, some false 
positives are expected. Follow-up of the candidates will screen out random fluctuations in 
the photometry. Identifying these transients as intrinsically high redshift supernovae will be 
facilitated by the deeper r and g-bands observations (24.7 and 25.0 magntiudes, respectively, 
for a single visit). 

If supernovae from massive Population III stars are still forming at 2; ~ 4, as suggested 
by Fig. m these sources will likely be well within LSST capabilities, thanks to a further 
reduction in the distance modulus and the increased sensitivity in the r band compared 
to the i band: These two factors will yield an overall gain of more than 1 mag, and a 5a 
detection of a SN2006gy-like source at z ~ 4 will require stacking of only 5 images in the 
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g-band. 

Our prediction for the observed pair instability supernova rate is plotted in Fig. [S] (for 
^wind = 60 km s~^) and shows that about one supernova per year per ~ 400 deg^ of sky is 
expected at 2; = 6, with the rate going down to one per ~ 1250 deg^ at z = 5. The rate at 
z ~ 5 is high enough to yield up to ~ 100 supernovae over the LSST lifetime. Extrapolating 
the rate for g-dropouts {z ~ 4) supernovae, we still expect a few events per year. 

Establishing the high-redshift nature of these transients is expected to be relatively sim- 
ple, as no other transient objects are known to possess such peak luminosities and light curve 
properties. A possible source of contamination arises fror n pair instability supernovae formed 



from the merger of two massive, metal-enriched stars (jPortegies Zwart fc van den Heuvel 



20071 ). Detailed estimates for the pair instability supernova rate from this channel are still 
lacking. Given the rarity of 2006gy-like supernovae in the local universe, the contamination 
is not expected to be severe, below 1 event yr~^. 

The detection of Population III supernovae will be more difficult if the progenitor stars 
are only moderately massive (typical masses 30 Mq). St ill up to 6 x 10^^ e rgs can by 
released by hypernovae with an initial mass of 20 - 25 M0 jNomoto et ahlboosh . However, 



the supernova rate will be higher, potentially tens of supernovae per deg^ assuming multiple 
stars per halo, and t he peak luminosity o f hypernovae might be fainter than that of a pair 



instability explosion (iNomoto et al.l l2003l ) . This will put them outside the reach of LSST, 
and a specific deep survey should be planned to search for them. First, a more detailed 
prediction for their light curve is needed. The best-case scenario for supernovae detection 
would be if multiple, very massive Population III stars formed per halo. 



Conclusion 



In this paper, we investigated the possibility of a long tail in the epoch of Population III 
star formation, extending to z < 7. This epoch has been suggested to extend to redshifts af- 
ter reionization because of inhomogeneous IGM metal enrichment and a finite transition from 
Pop I II to Pop II modes of star formation fe.g.. IScannapieco et al.ll2003l : iJimenez fc Haiman 
20061 : iTornatore et al.l 120071 : IShuU fc Venkatesanll2008l ). Our simulations reach a mass reso- 
lution sufficient to identify Population III star formation driven by H2 cooling at z < 26 in 
halos with M > 10^ Mq. At the same time, they cover a sufficiently large volume (10^ Mpc^) 
to be representati ve of the Universe a t z > 5 for the halo mass scale relevant to Population 
III formation (see iBada fc Ravll2005h . Coupling the dark matter only dynamics with con- 



servative assumptions on self-enriched (genetic) and wind enriched scenarios, we show that 
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pockets of metal-free stars clearly exist well after the end of reionization, although with a 
small number density (one new halo Gpc~^ yr~^ formed a.t z = Q ) . Ou r results confirm the 
Population III star formation trend identified by lTornatore et al.l (120071 ) for z <6, although 
their mass resolution was insufficient to evaluate self-enrichment from minihalos, which might 
explain the different behaviour at 2; > 6 (see Appendix However, our inferred Population 
III star formation rate assumes one very massive st ar per halo and is sma ller by a factor 
~ 100 at 2; ~ 6 compared to that plotted in Fig. 1 of iTornatore et al.l (120071 ). If we were to 
assume instead that a fraction e ~ 10^^ of the metal-free gas in converted into P opulation 
III stars, then our results would be much closer to those of ITornatore et al.l (120071 ) at z < 6. 



While the 2; < 10 Population III star formation rate is several orders of magnitude lower 
than that of metal-enriched stars, the presence of metal-free gas at z < 6 opens an interesting 
window for detecting Population III stars. If their expected initial mass function is peaked 
toward very massive stars (M ^ 100 Mq) that could thus explode as hypernovae or pair- 
instability supernovae, their expected rate (see Fig. [5]) and peak luminosity of ttiab ^ 26 
at 2 < 5 suggest that future ground-based searches with LSST can probe this mode of 
star formation. We predict one z ^ 5 Population III supernova per year per 1250 deg^. We 
suggest that large-area blind searches will be necessary to find these objects, as the formation 
sites for Population III stars are distributed almost uniformly on the sky, far (~ 150 kpc 
equivalent to ~ 5 arcsec) from the nearest (dwarf) galaxy. A prediction of our model for 
metal-free supernovae is that no host galaxy will be present. Furthermore, the two point 
correlation function of these supernovae will show that these sources are significantly less 
biased than the general population of ~ IO^Mq dark matter halos at that redshift (see 
Fig. [3D. 

All our assumptions about the enrichment of the IGM and the formation rate of metal- 
free halos are robust. Thus, the presence of pristine gas available for star formation at 2; < 6 
is a very strong prediction of our modeling. However, translating the amount of gas available 
into an expected supernova rate carries a major uncertainty, namely the Population HI initial 
mass function and the star formation efficiency. These are highly debated topic, in particular 
the IMF of Population HI stars. In fact some tension betw een theory /si mulations and 



observ ations is present (iTumlinson et al.ll2004j : lTumlinsonll2006l ). In addition. lYoshida et al. 



(I2OO7I ) found with a 1-D model for cooling of metal free gas in the presence of a strong UV 
background, that the IMF of Population HI stars is shifted toward smaller masses compared 
to formation at higher redshift. Under such conditions, it is unclear whether fragmentation 
will lead to multiple Population HI formation. Properly addressing the cooling process of 
these halos is therefore a fundamental step toward evaluating the Population HI formation 
rate at the edge of the reionization epoch. These issues will be examined in a follow-up 
paper. 
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A. Validation and Convergence 

To correctly evaluate the number density of metal-free halos around the epoch of reion- 
ization, it is necessary to ensure that the numerical simulations we used have a sufficient 
resolution. Strong winds can originate only from relatively massive halos, and are thus 
relatively easy to identify. However, self-enrichment is more difficult to model because it 
may arise primarily from minihalo progenitors. Identifying these halos requires a high mass 
resolution (single particle mass below ~ 5 x 10^ Mq). 

A strong hint that our numerical results have reached convergence with respect to self- 
enrichment comes from the excellent agreement of t he Population III h a lo for mation rate with 



the analytic model for self-enrichment presented in lTrenti fc Stiavellil (120091 ) (see Fig. [2]). In 
addition, we compare here our production run with N = 1024^ and a box size of 7 h'^ Mpc 
(single particle mass of 3.4 x 10^ M0) to the results of a control run at lower resolution, 
with N = 512^ and box size of 10 Mpc (single particle mass of 7.9 x 10^ Mq). In Fig. [H] 
we plot the corresponding halo formation rate for our reference model with a wind speed of 
60 km s~^. We clearly see that at 2; < 7 the two simulations yield the same results. The 
relative difference in d?n/dtdV is below the variance of the halo number counts in a single 
snapshot. Note that at higher redshift (2; ~ 15) the lower resolution run exhibits a peak in 
Population III formation rate, which is a numerical artifact induced by the inability to follow 
self-enrichment due to minihalos formed at z > 20. The resolution of this "low resolution" 
control run is equal to the highest resolution run of Tornatore et al. (j2007 ) suggesting that 
their results might not be reliable at 2; > 10. 
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Fig. 1. — Black solid line: minimum halo mass required for cooling and collapse of a metal- 
free halo in presence of an H2 dissociating Lyma n - Wern er background. Calculations are 
based on the standard model by lTrenti fc Stiavellil (120091 ) — see their Fig. 1. For z < 13.4, 
molecular hydrogen cooling is strongly suppressed in minihalos, and Population III star 
formation occurs in halos with virial temperature Tvir = 10'^ K (dashed blue line). At higher 
redshift, cooling is instead possible in minihalos with a minimum mass MjninH2 (^^d dotted 
line) self-consiste ntly derived from the Ly man- Werner background created by the first stars 
and galaxies (see lTrenti fc Stiavellil l2009l for details of the model). 
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Fig. 2. — Solid black line: formation rate of halos that have not experienced self-enrichment 
(number of new halos formed per proper unit time (yr) per unit volume (comoving Mpc^)). 
The red dashed line is the pr ediction for the sanie qu antity based on the Extended Press- 
Schechter ("EPS") model of iTrenti fc Stiavellil (120091 ) . There is a remarkable agreement 
especially at redshift z < 10. Winds are not taken into account in this figure. 
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Fig. 3. — Comoving distances and masses of the nearest-neighbor halo of non self-enriched 
(red dots) and truly pristine halos (black squares) at redshift z = 6 for a wind speed fwind = 
60 km s~^. A significant fraction of halos that were not self-enriched have neighbors too 
distant to pollute them with their metal winds. The blue dashed line gives the average 
distance of halos of mass M under the assumption of a uniform distribution, highlighting 
the effect of clustering. Truly pristine halos are significantly less clustered than their metal 
enriched counterparts. 
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Fig. 4. — Population III halo formation rate after accounting for wind pollution (short 
dashed green line fwind = 40 km s"^; solid black line f„ind = 60 km s~^; long dashed blue 
line fwind = 100 km s~^) compared to the formation rate of halos that were not self-enriched 
(dotted red line). The red dotted line is also equivalent to assuming fwind = 0. 
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Fig. 5. — Observed Population III supernova rate per year per deg^ as a function of redshift 
for our standard model of wind propagation (wwind = 60 km s~^). This has been obtained 
assuming one Population III supernova per halo. The field of view depth is = 1 centered 
at redshift a choice representative for typical Lyman-Break dropout searches. 



-24- 




Fig. 6. — Population III halo formation rate obtained with our standard wind prescription 
("i^wind = 60 km s~^) from two cosmological simulation, a "Low-Resolution" = 512^ run 
(solid blue line) and a "High- Resolution" = 1024'^ run (solid black line). The dotted blue 
and black lines represent the rate in presence of self-enrichment only. The results of the two 
runs agree well for < 8. At higher redshift (8 < -z < 17), the low resolution run is not able 
to follow self-enrichment and thus over-estimates the formation rate of metal-free halos. 
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Table 1: Sensitivity of current and future facilities for Population III detection 



Telescope /Instrument 


Band 


Sensitivity 


Integration Time 


Field of view 


Time/ (10^ deg^) 


LSST 


i 


rriAB = 26 


1200 s 


9.6 deg^ 


35 h 


SUBARU/Suprime-Cam 


i 


ruAB = 26 


470 s 


0.25 deg^ 


522 h 


HST/ACS 


F775W 


rriAB = 26 


800s 


11.3 arcmin^ 


7 X 10^ h 


JWST/NIRCAM 


F090W 


niAB = 26 


< 30 s 


9.7 arcmin^ 


3093 h 


TMT/IRIS" 


i 


niAB = 26 


~ 30 s 


0.06 arcmin^ 


5 X 10^ h 


JWST/NIRCAM 


F090W 


niAB = 31 


65 h 


9.7 arcmin^ 


N/A 


TMT/IRIS" 


z 


rriAB = 31 


34 h 


0.06 arcmin^ 


N/A 



Note. — Summary of the sensitivity for 4.5 ^ z < 7 Population III SNe searches using the current 
and future facihties listed in the first column. The second column reports the filter used, the third the 
5cr sensitivity reached in the time given in the fourth column. The fifth column gives the field of view 
and the last column the time needed to survey 10'^ deg^: SNe detection is maximally efficient with LSST 
given its superior field of view. Last two lines: a deep field with JWST NIRCAM or with a 30 m class 
telescope will approach the depth necessary for direct imaging of a cluster of Population III stars, if these 
exist, a: TMT sensitivity is an estimate based on the "TMT Instrumentation and Performance Hand- 
book" fhttp: / /ww w. physics. uci.edu/TMT-Workshop/TMT-Handbook.pdf| ), as no exposure time calculator 
is currently available. 



